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We study boson star configurations with generic, but not non-topological, self-interaction terms, 
i.e. we do not restrict ourselves just to consider the standard A|'i/'|* interaction but more general 
U(l)-symmetry-preserving profiles. We find that when compared with the usual potential, similar 
results for masses and number of particles appear. However, changes are of order of few percent 
of the star masses. We explore the stability properties of the configurations, that we analyze using 
catastrophe theory. We also study possible observational outputs: gravitational redshifts, rotation 
curves of accreted particles, and lensing phenomena, and compare with the usual case. 



I. INTRODUCTION 

In the past two decades, the study of the universe in its 
earlier stages has given an important role to scalar fields: 
conventional models for inflation rely on its possible po- 
tential energy and topological defects may form when a 
scalar field breaks some fundamental symmetry. Both of 
these processes may be responsible for the formation of 
the large scale structure we now see. In addition, scalar 
fields connected with time or space-time variation of fun- 
damental constants have produced the most powerful al- 
ternative theories of gravitation, known as scalar-tensor 
theories, and entered in all unification scheme trials made 
so far. It must be said, however, that no fundamental 
cosmological scalar field has been directly observed yet. 

If one accepts the need of scalar fields in the cosmo- 
logical scenario, one interesting question naturally arise: 
may these scalars be the seed of astrophysical structures 
or of observable phenomena that could signal their exis- 
tence? Objects made up of scalar massive particles were 
introduced as early as 1968 by Kaup and Ruffini and 
Bonazzola These configurations are now known as 
boson stars. These stars, contrary to the more common 
neutron or fermion stars, are not supported by Pauli's ex- 
clusion but by Heisenberg's uncertainty principle, which 
effectively keeps scalars from being localized to within 
their Compton wavelength and prevents their collapse to 
a black hole. 

The first investigations with a nonlinear "0^ potential 
were carried out by Mielke and Scherzer They cal- 
culated this potential from a nonlinear Heisenberg-Pauli- 
Weyl spinor equation and found the first scalar field so- 
lutions with nodes. More recently, Colpi et al. Q proved 
that the existence of a self-interaction in the boson La- 
grangian could yield higher values for the masses of the 
configurations. This means that for some values of the 
free parameters existing in the theory, stellar structures 
of appreciable masses and extremely high density may 
arise. This could have incredible effects on our under- 
standing of the non-baryonic content of the universe and 
have different observational effects: affecting usual ob- 
jects as galaxies j|] or even light trajectories in powerful. 



degree-ranged, microlensing phenomena [|6[. These kind 
of features, as well as the important point of star rota- 
tion (see 0] among others), was thoroughly studied in 
the past few years, including gauge charged stars, boson- 
fermion models (see 1^ for reviews) and models of stars 
in alternative theories of gravity . 

In all these works, the self-interaction was always 
choiced to be like Aj-i/'l^, thus giving a matter sector pro- 
vided by, 

An = -^S^'' d^rdA - im^lV'P - Ja|V|4 , (1) 

where A is a constant. However, a priori, there is no 
reason to maintain just the form adopted for the self- 
interaction and not a more generic potential term T^(V'). 
Although in order to retain the U{1) symmetry of the 
whole Lagrangian, this potential must be a function of 
It/'P, in principle, it may have any other particular form. 
In this work we would like to explore whether substan- 
tial variations in the form of the potential may yield ap- 
preciable changes in the final configurations, the binding 
energies, or the masses of boson stars. We also analyze 
the changes that a different potential introduces in grav- 
itational redshifts and microlensing phenomena. This 
could be important to discover if different scenarios for 
dark matter candidates, explanation of galactic rotation 
curves H, o r new observational gravitational effects may 
appear [|0| . Additionally, new potentials may strongly 
affect gravitational memory and evolution scenarios pd| ] 
within theories with time variation of Newton's constant 
(for the most recent account of this see Ref. [p^). 

Knowledge about that modifications in the boson po- 
tential may yield to profound changes in the structure 
of the stellar objects comes from the works of Friedberg, 
Lee, and Pang |]lj-|l^ on non-topological soliton stars. 
However, we point out that a systematic study on their 
possible observational signatures is still absent. We shall 
briefly mention their features and compare with our po- 
tential choices below. However, our main aim here is to 
study ij potentials which are not non-topological ones may 
still introduce significant changes in the configurations. 

This work is organized as follows. In the next section 
we introduce the models for boson stars and search, with 
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different self-interaction terms, for their configurations 
and stability properties. Section III compares the results 
of these generic boson star models, concerning observa- 
tional outputs, with those obtained for the usual A|i/)|^ 
term. We give our conclusions and provide a brief dis- 
cussion in the final section. 



II. THE MODELS 

We shall take into account the following boson La- 
grangian 



2" ^ 
with the potential 



Ui\^P\^)=m^\^\^ + -V{\M^) 



(2) 



(3) 



This action possesses an invariance under the global C/(l) 
transformation tp e^^ip, which gives rise to a conserved 
current 

= tg'^^' - i'd^r) , (4) 

and a corresponding conserved charge 

N = J dx^y/^J°. (5) 

This conserved quantity can be identified with the num- 
ber of bosons present in the stellar structure. We shall 
consider the usual spherically symmetric metric 



-B{r)dt^ + A{r)dr^ + r^dfl^, 



(6) 



and shall also demand a spherically symmetric form for 
the field which describe the boson, i.e. we adopt the 
ansatz: 
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and the Klein-Gordon equation for -0, 
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In these equations, we have used dimensionless units, 
which are common to those introduced in the paper by 
Colpi et al., 



X = mr , 



(11) 



for the radial coordinate ( ' stands for the derivatives with 
respect to x) and. 



PI 



(12) 



where Mpi = G^^/^ jg the Planck mass. 

In order to consider the total amount of mass contained 
within a radius x we have changed the function A in the 
metric to its Schwarzschild form, 

Ai.)^[,-'jfiy\ ,13, 

Note that the terms corresponding to the potential are 
correctly normalized. From the Lagrangian (^) one may 
see that V has dimensions of [Energy]^, and all of them 
are divided by two squares of masses. However, in order 
to avoid the explicit appearance of the boson mass m and 
the Planck mass Mp\, we need to define the form of V. 
We shall look at several choices. 



ip{r, t) — xif ) exp [— itui]. 



(7) 



This ensures that a static solution is still possible. To 
establish that this time dependence is the lowest energy 
solution at a fixed number of particles, it is necessary to 
make a first order variation 5{E — wN), where in general 
is a Lagrangian multiplier associated with the conser- 
vation of N and E is the energy (mass) of the system. 
See for instance the appendix of Ref. ||lj] for a detailed 
derivation. This form for ■0 is independent of the poten- 
tial U . For a more detailed derivation of the basic model 
we refer the reader to Refs. Using the Einstein 

gravitational action we obtain the field equations. 



2M' 



B 

1 / 47r 
2 



A 



y(crMpi/V47r), 



(8) 



The first expansion to nonlinear potentials can be 
found in the 1981 paper by Mielke & Scherzer They 
constructed a potential for the Klein-Gordon equation 
from the Heisenberg-Pauli-Weyl nonlinear spinor equa- 
tion. It has the general form U = m^j^p _ q,j^|7^|4 _|_ 
a2|'0|^, where ai and a2 are two positive constants. They 
presented solutions with nodes for the first time. 

The standard (CSW g) choice is V{ij) = A|^/'|^ with A 
a constant. Here, the usual adimensionalization appears, 



47r 



477 



t/(CTA/pi/V47r) 



(4^) 
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(14) 



where A = AA/^i/4 nm . As we stated in the intro- 
duction, with this choice, the order of magnitude of 
boson star masses is deeply enhanced. It grows from 
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M - MlJm when A = to M MIj/to^ when A 7^ 0. 
Recah that the mass of a neutron star is roughly given by 
the Chandrasekhar mass Mch ~ -^pi/'^n which is close 
to a solar mass (to„ is the neutron mass). 

Other three options we would like to explore are (note 
that these are options for t/d?/;!^), given in Eq. (j^), not 
only for as will be clear below): 

• Cosh-Gordon potential: 



Then, in order to be consistent with equations (^) and 
(p|-|l0|), the field equations and the definition of U , and 
to avoid a useless double counting of the mass term, we 
must take particular choices for a; the parameter (3 is 
still free for choice. 

• Cosh-Gordon potential: 

a ^ 2(M|i/47r)//32 



C/cosh = am^ [cosh(/3/^) - l] 
• Sine-Gordon potential: 

C/,in = am^ [sin(^/2[/3/^ - 1]) + l] 



Sine-Gordon potential: 



'[1 



cos 



• t/(l)-Liouville potential: 

C/e,p = am2[exp(/32|V;n-l] 

The usual Liouville potential exp{(3ip) has to be 
changed so that a t/(l) symmetry is ensured. 

Let us first consider a series expansion of these poten- 
tials. In order to do so we shall consider a value of f3 
such that, when going from -0 to the dimensionless cr, 
the arguments of the functions are not affected. The pa- 
rameter P is arbitrary, it enlarges the parameter space 
of the solutions, as was the case with A in CSW's so- 
lutions. The appearence of the factor f3 is just because 
dimensional grounds, while a can be used to get a sim- 
ple mass term in the boson Lagrangian, as we shall see 
below. Taking this into account, the series expansions 
are. 



• J7(l)-Liouville potential: 

a = {Mpi/Att)/P^ 

In this way, the potential V is everything but the first 
terms in each of the previous series. Hence, it is a series 
of attractive-repulsive self-interactions in the case of the 
Sine-Gordon potential and a series of repulsive poten- 
tials in the Cosh-Gordon case. The case of the U{1)- 
Liouville potential is reminiscent of the Cosh-Gordon 
one, in the sense of being a series of repulsive power law 
self- interactions, just the coefficients differ. The form of 
these potentials and others mentioned above are shown 
in Fig. 1. 

For the numerical procedure, it is best to make a re- 
definition of the scalar field mass 



= amn? 



(18) 



and with this also redefine the coordinate x and the fre- 
quency ri. Notice that /? still appears within the differ- 
ential equations while a does not. 



Uc.„f.h = am 
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= am X 
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(15) 



(16) 



(17) 



Note that, from each expansion, we are recognizing a 
usual mass term (proportional just to m^). This term 
is very important: without it, it is impossible to find 
solutions with exponential decrease of the scalar field, 
something relevant for the definition of the star radius. 



A. Soliton stars 

A potential with symmetry breaking was investigated 
by Lee et al. |p^-|l5|. They called the solutions non- 
topological soliton stars, and found that the mass has 
units of Mp[/(mag) which is huge in comparison with a 
boson or neutron star (for the case of comparable boson 
and fermion masses). The potential investigated was U = 
TO^|V'P(1 — \4'\'^/'^o)^ where ctq is a constant; this belongs 
to the more general forms of potentials derived in 

Compared with the usual boson star case, non- 
topological soliton stars have to fulfill two characteristics: 

1. The Lagrangian must be invariant under a global 
U{1) transformation. 

2. In the absence of gravity, the theory must have 
non-topological solutions; i.e. solutions with a fi- 
nite mass, confined to a finite region of space, and 
non-dispersive. 
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In general, boson stars accomplish the requirement 1. but 
not 2. Invariance under U{1) only requires that the po- 
tential be a function of and in order to fulfill con- 
dition 2., U must contain attractive terms. This is why 
the coefficient of of Lee's potential has a nega- 

tive sign. Finally, when 1-01 — + oo, U must be positive, 
which leads, minimally, to a sixth order function of ip for 
the self-interaction. It is then clear that CSW's, J/cosh, 
and C/exp choices are not non-topological potentials. Nei- 
ther of them have attractive terms. This is why, a priori, 
we may say that the order of magnitudes for the boson 
star masses remains the same as in CSW's case. The 
Sine-Gordon potential Usin has, on the contrary, a simi- 
lar series expansion, up to the sixth order, to that corre- 
sponding to Lee's potential. But here, what happens is 
that [/sin, in the absence of gravity and for a real scalar 
field, has not a non-topological soliton solution, instead, 
it has a topological one. It has a degenerate vacuum: an 
infinite set of a values for which {/gin = 0. For a detailed 
account of this, we refer the reader to Lee's book, espe- 
cially Chapter 7 and exercise 7. 1 . Then, [/gin neither 
is a non-topological soliton potential. 

B. Numerical solutions 



Ms = 47ry pr^dr (19) 



^ ^ 

A second mass formula can be derived for a general quasi- 
static isolated mass insula where a time-like Killing vec- 
tor — 2n" exists 0. Tolman's expression is: 




For an asymptotically flat spherically symmetric space- 
time, both masses agree with each other. 

Finally, Figs. 6 and 7 shows the behavior of the U{1)- 
Liouville potential. It is both, qualitatively and quanti- 
tatively similar to the usual CSW's case. We have stable 
and unstable branches. The maximal mass and particle 
number are higher by about 5% which are the largest de- 
viations from the standard case. The repulsive potential 
terms yield larger contributions in comparison with the 
Cosh-Gordon potential. 



Fig. 2 shows the usual plot of boson star configurations 
for the case of the Cosh-Gordon potential. The maximal 
mass is slightly higher than in the standard case due to 
the additional higher-order repulsive terms in the poten- 
tial. For /3 = 1, we find M„iax = 0.638 M^i/m and 
N,na^ = 0.658 MlJrn^ what for M and N is higher by 
0.5%. Fig. 3 shows the stability analysis, which can be 
done using catastrophe theory |17|-p^. One necessary 
condition for the configurations to be stable is a negative 
binding energy. However, this is not sufficient. Fig. 3 
shows the appearence of two cusps signaling that there is 
a change in the star stability. The first branch is the only 
stable one, while the second and third are both unstable. 

Fig. 4 represents similar profiles, but for the Sine- 
Gordon potential. In this case, the maximal values of 
mass and particle number are below the pure mass po- 
tential case. The influence of the higher order attractive 
terms is noticeably. For M and iV, the maximal val- 
ues are lower by about 2%: Afmax — 0.620 M'^^/m and 
-^max = 0.639 Mp^/m"^. Fig. 5 shows the bifurcation 
plot for this case. The diagram in Fig. 5 shows cusps 
again, where is a change in stability. A remark on the 
calculation of mass is in order. We check that the cal- 
culation of the mass is correct by applying two different 
mass definitions. The first is the Schwarzschild mass, 
which is defined by the energy density p and which also 
appears in the asymptotic spherically symmetric space- 
time, B{r) ^ 1 — 2M/r, where M is the mass of the 
boson star. The formula for the Schwarzschild mass is 



III. OBSERVATIONAL OUTPUTS 

At this stage, we have succeeded in proving that a 
change in the form of the self-interaction among the 
bosons yields appreciable -but small- changes in the form 
of the star configurations. We shall discuss below the fea- 
sibility of detecting observational consequences of these 
results concerning gravitational phenomena. In particu- 
lar, we are interested in seeing if appreciable differences 
appear in the computation of gravitational redshifts, ro- 
tation curves, and gravitational lensing features. For the 
Sine-Gordon potential, and because of the smaller masses 
it produces, gravitational phenomena are diminished. So 
we shall be mainly interested in Cosh-Gordon and U{1)- 
Liouville cases. 



A. Gravitational redshifts 

In this section, we follow Ref. jl^], and make use of the 
assumption that the scalar particles have no interaction 
-other than gravitational- with baryonic matter. Thus, 
this normal matter can penetrate the boson star up to 
the center and if it emits radiation there, well within the 
gravitational potential, we expect the spectral features 
to be redshifted. 

The gravitational redshift z of our static line element 
is given by 

1 + ,= ,S4, (22) 
Y B{int) ' ^ ' 
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where int stands for the position of the emitter particle 
with respect to the star center. As the receiver is prac- 
tically at infinity, B(oo) ~ 1. The maximum possible 
redshift is obtained when the particle emits exactly at 
the center of the boson star, where the metric deviates 
maximal from outside vacuum space-time. 

We are only interested in stable configurations, the 
maximum redshift is then provided by the maximum 
value of cr(0), which gives the biggest mass. As it was 
shown in Ref. the simple mass term produced a 

maximum redshift of 0.45 while for CSW's choice, with 
A tending to infinity, one gets 0.69. Here we find Zmax = 
0.46 for the Cosh-Gordon potential and z^ax = 0.49 for 
the t/(l)-Liouville potential. For comparison, we quote 
the result z^ax = 0.47 for neutron stars. 



B. Rotation curves 

Another gravitational effect considered for CSW's 
choice [0, and which we would like to compare with 
the more generic potentials here studied are the rotation 
curves of test particles moving around boson stars. For 
our metric, circular geodesies have a rotation velocity (as 
measured by an observer at infinity) given by. 



xB' 



(23) 



In Fig. 8 we compare the rotation curves for the case 
A = 0, A = 100 of CSW's choice, and our potentials: 
Cosh-Gordon and [/(l)-Liouville. Already from the usual 
case, it was shown that the possible velocities that parti- 
cles can reach are a notable amount of the speed of light, 
and matter can have an impressive kinetic energy. If such 
kinetic energy were transferred to radiation, we could 
expect very luminous boson stars, orders of magnitude 
more luminous than the Sun. This provides speculative 
alternatives to accretion disks around black holes, and 
would make boson stars almost indistinguishable from 
its final effects. This is currently being analyzed, espe- 
cially concerning the possibility of having boson stars at 
the center of some galaxies, as was proposed with neu- 
trino balls. For the new potentials we are analyzing here, 
we obtain that C/cosh and C/oxp produce similar angular 
velocities to the A = case. Their maximum velocity 
happens for values x ~ 5 and have typical magnitudes of 
100 000 km s-^ 



C. Gravitational lensing 

Boson star microlensing effects were first investigated 
by D§,browski and Schunck ||], for CSW's A = case, 
also known as the mini-boson star. Their procedure, 
which we closely follow, consist in studying the pho- 
ton trajectories along the curved (boson star generated) 
space-time. For particular details of the derivation of 



quoted formulae see Refs. §,|l|. In Ref. a related 
study on microlensing features was made, taking a scalar- 
field-generated naked singularity as lens. It has the prop- 
erty of producing both, positive and negative binding an- 
gles; in this later case, in a way similar to the recently 
studied wormhole microlensing scenario . 

The light traveling from a distant source is deflected, 
because of the presence of the boson star, with a deflec- 
tion angle given by (see Fig. |^ for a schematic drawing 
of the geometry). 



a = Aip — TT, 



where, 



dr 



(24) 



(25) 



with impact parameter b — tq^JX/ B{rQ) (rg is the clos- 
est distance between the light ray and the center of the 
boson star: the first point where the square root in the 
denominator is non- negative) . The lens equation can be 
expressed as the difference between the true angular po- 
sition, and the image positions, "iJ, as P3,M 



sin(?9 — (3) = sind 



(26) 



where Dps (Dos) stands for the angular distance between 
the point P close to the lens and the source (the observer 
and the source). Also, from the geometry of the lens we 
have sin?? = b/Doi- Hence, choosing ■& and the distance, 
we have b, and Aip may be computed afterwards. In the 
numerical program, we use again the dimensionless quan- 
tities of Eqs. (11-12) and instead of the impact parameter 



6, we follow M and take The term {r/bY in (25) is 



then x^/{w'^Sli sin^ -d) or just x'^ / {vj"^ Dlid"^) for small i?, 
respectively. Our numerical program uses always the cor- 
rect sin I? without any abbreviations so that also angles 
in the degree regime can be calculated. The examples in 
Figs. ^ and ^ apply d in arc-seconds having the addi- 
tional unit factor 1/206265 for one arc-second in radians. 
The change to other units can then be explained by an 
additional distance factor n. For instance, if d has to 
be measured in milli-arc-seconds, n equals 10^^. In or- 
der to get rid of the distance factor within the numerical 
program, it is chosen Doi — 206265/(tx7n). 

Furthermore, the reduced angular deflection angle is 
defined to be 



a — "d — B — arcsin ( — — sin a 

\Dos 



(27) 



A second lens equation can be derived for large de- 
flection angles, where Dis cannot be considered as being 
similar to Dpg. Of course, by its construction the source 
is always within a plane with constant distance to the 
observer, and studying the diagramatic view depicted in 
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Fig. 9, it can be obtained that (for a detailed derivation 
see Appendix of Ref. M), 



Dis , 
sm a — — — cos v cos 



arcsin ( —2^ sin(i? — a) 



X [tani? + tan(Q! — -d)] , 



(28) 



where Dis stands for the angular distance between the 
lens and the source. Lens equation ( |2^ ) requires only the 
proportion between Dps and D^s so that, in general, the 
position of source S describes a more complicated surface. 
For the physical situation considered in our paper, the 
differences for a amounts a few parts per thousand of a 
degree at most so that our Figs. 10 and ^ describe both 
cases. 

Assuming that the boson star lens is half-way between 
the observer and the source, such that Dis/Dos = 1/2 
and Dps /Dos = 1/2, we performed numerical computa- 
tions of the reduced deflection angle for our new poten- 
tials, which we show in Figs. |l^ and [l^. The difference 
among these cases and the simple mass term (correspond- 
ing to the mini-boson star) is clearly observable. We have 
taken for the plot the maximum central density (which 
produces the maximum deflection angle). In the case 
of the mini-boson star, the biggest possible value of a 
is 23.03 degrees with an image at about 1) = n x 2.88 
arc-sees with the distance factor n — n{Doi,^l) which is 
a function of the distance from the observer to the lens 
and the scalar field frequency, whose inverse can be asso- 
ciated with the star radius. In our examples we assume 
that n = 1, which fixes "d to be measured in arc-sec. 

The characteristics of the boson stars produce the 
deflection angles which depend on the observer-to-lens- 
distance Doi and the mass. As mentioned above, if 
is chosen to be of the order of arc-sees (distance factor 
n = 1), then the distance Dqi is measured in units of 
206265/ci7. Under the assumption that the mass of the 
boson star is lO^^M© one has m ~ 10~^^cm~^. Then, 
the distance Doi is about lOOpc. If the distance factor 
is n = 10"'^, and so measured in milli-arc-secs, the 
boson-star-lens is at about lOOkpc. 

We assumed that the boson star is transparent and cal- 
culated the deflection angles. All qualitative features of a 
non-singular spherically symmetric transparent lens can 
be revealed using Figs. ^ and |l^: the lens curve for the 
maximal boson star. Three images exist, two of them 
being inside the Einstein radius and one outside. An 
Einstein ring with infinite tangential magnification (tan- 
gential critical curve) is found, and also a radial critical 
curve for which two internal images merge. The appear- 
ance of the radial critical curve distinguish boson stars 
from other extended and non-transparent lenses. For a 
black hole or a neutron star, the radial critical curve does 
not exist because it is inside the event horizon or the star. 
Two bright images near the center of the boson star and 
the third image at some very large distance from the cen- 
ter are found. For non-relativistic cases smaller angles 
will be found. An interesting point can be made if one 



considers an extended source. In such a case one finds 
the two radially and tangentially elongated images very 
close to each other. Then, looking along the line defined 
by these two images the third one can be detected at a 
very large distance. 

For the Cosh-Gordon potential, we obtain as maximal 
reduced deflection angle 23.229 degrees, and for C/cxp, an 
even larger deviation, the biggest possible being 24.391 
degrees with an image about at the same place. Differ- 
ences among these cases and the usual one is between 0.2 
and 0.4 degrees. 



IV. DISCUSSION 

In a recent communication. Ho et al. have studied 
the maximum masses of boson stars formed with different 
self-interaction terms (all of them, however, of power law 
form and of positive sign). By just comparing the order 
of magnitude of the terms involved in the self-interaction 
with the mass term, and asking for them to be of the 
same order, they were able to see that the contribution of 
higher order terms goes as power of (m/Mpi)^. The fact 
that the contribution to the maximum masses of these 
different boson stars decreases (if m < Mp\), does not 
automatically yield to unobservable effects, as we have 
discussed in the previous sections. 

The star masses maintain the order of magnitude, for 
equal single boson masses, when compared to those cases 
studied by CSW B, which is in agreement with the re- 
sults of Ho et al. [Ea. This also stems from the fact that 
all Lagrangians analyzed are not non-topological ones. 
Changes are of order of several percent of the star mass. 

For the first time, we have investigated the systems of 
differential equations of Einstein-Cosh-Gordon, Einstein- 
Sine-Gordon, and Einstein-J7(l)-Liouvillc. The different 
potentials studied so far showed similar gravitational red- 
shifts and rotational curves, with high angular veloci- 
ties, when compared among them, and between them 
and the CSW's case. However, large deviations in the 
maximum angular deflection have appeared, with differ- 
ences amounting appreciable parts of a degree. Some 
observational effects distinguish boson stars from other 
non-transparent compact objects. 

However, fair is to say that if an observable determi- 
nation proves the existence of a boson star, the effec- 
tive form of the Lagrangian may be hidden within the 
percentage of possible errors. In that case, facing with 
the problem of degeneracy -i.e. different physical theo- 
ries giving the same observational effects- Occam's razor 
would probably lead us to consider just the CSW's choice. 
Only a detailed knowledge of the boson involved, and the 
average form of the interactions, all of them encompassed 
in the self-interaction term, may shed light on the explicit 
model for the Lagrangian. 

The definition of the actual boson which participates in 
the construction of the boson star may also influence the 
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transparent consideration for gravitational lensing phe- 
nomena. For instance, if the star is made up of Higgs 
particles, we may expect some kind of interaction apart 
from the gravitational one that may yield the star to be 
a non-transparent object. 

Furthermore, we become aware of the possibility that 
some potentials may give rise of tunneling of parts of the 
scalar field. Of course, this effect can only occur in the 
quantum regime, hence, if boson stars are in the order 
of magnitude of atoms or even atomic nuclei. Our first 
preliminary results for the Newtonian case show that es- 
pecially the form of potentials of Lee et al. and of Sine- 
Gordon can lead to instability due to tunneling. The 
effect could mean two things: (i) The boson (soliton) 
star is destroyed: it disperses or it forms a black hole, 
(ii) The boson (soliton) star experiences an internal re- 
arrangement. We expect to report on these issues on a 
forthcoming article. 
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KRB 
Cosh-Gordon 



CSW 
Sine-Gordon 
MS 



U(1)-Liouville 

FLP 

FIG. 1. Comparison of different potentials. KRB 
(Kaup-RufRni-Bonazzola): , Cosh-Gordon: 2(cosh(a;) — 1), 
?7(l)-Liouville: exp(a;^) - 1, FLP (Friedberg-Lee-Pang): 
x^il - x'^/2f, CSW (Colpi-Shapiro-Wasserman): 
+ imx*/2, Sine-Gordon: (8/7r2){sin(7r/2[a; - 1]) + 1}, MS 
(Mielke-Scherzer): x^ - x"^ + x^ /2,. 
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FIG. 2. Configurations for I5?^&i stars self-interacting via 
a Cosh-Gordon potential for /3 = 1. 
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FIG. 3. Stability analysis for Cosh-Gordon configurations. 
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FIG. 4. Configurations for boson stars self-interacting via 
a Sine-Gordon potential. 
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FIG. 7. Stability analysis for the L''(l)-Liouville poten- 
tial configurations. The binding energy changes sign and the 
presence of cusps signals changes in stability. The first branch 
is the only stable one. 
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FIG. 8. Rotation curves for non-interacting boson stars 
and boson stars with generic potentials. Shown in the Figure 
are curves corresponding to A = 0, 300 of the CSW's choice, 
taken from the work of Schunck and Liddle, and our new 
results for the potentials Ucoah and Ucxp- The maximal veloc- 
ities are: 122 990km/s at a: = 20.1 for A = 300, 102 073km/s 
at X = 4.1 for A = 0, 104 685km/s at x = 4.2 for Ue^p, and 
102 459km/s at a; = 5.9 for Ucosh- 
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FIG. 9. Schematic lens diagram for microlensing phe- 
nomenon. Dol is the distance from the observer (O) to the 

lens (L), Dos is the observer-source distance, and Dis the dis- 
tance from the lens to the source (S). The angle j3 denotes 
the true angular position of the source whereas & measures 
the angle of the image position. Dps is the distance from the 
point P in the source plane to the source (for small deflec- 
tion angles, Dps ~ Dis, as well as Dos ~ Doi + Dis). For 
the meaning and use of all other parameters see text and the 
Appendix of the paper by D^browski and Schunck. 



FIG. 10. The reduced deflection angle (difference between 
the true and the image angular position) as a function of 
the the image position for the different potentials, f/exp pro- 
duces the largest deflection angle in comparison with the sim- 
ple mass term potential (CSW's choice with A = 0) and the 
Cosh-Gordon potential. We have chosen here the stable max- 
imal mass boson star. 
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FIG. 11. Comparison for all three potentials at the maxi- 
mal reduced deflection angle. The upper curve is for J7exp, the 
middle curve is for J/cosh and the lower one for A = CSW's 
choice. 
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